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Abstract 


The  nonlinear  transmission  of  Znrtetrabenzporphyrin  was  measured  in  a 
Z-scan  setup  with  532-nm  wavelength  laser  light  and  a  13-ns  pulse  duration. 
The  excited-state  absorption  cross  section,  the  excited-state  refractive  index 
cross  section,  and  the  linear  and  nonlinear  absorption  contribution  to  a 
thermal  index  change  were  investigated.  The  effects  of  fluorescence  and 
acoustic  waves  on  the  nonlinear  response  of  TBP  have  been  determined. 
Limiter  performance  was  modeled  in  an// 14  limiter,  and  saturation  effects 
were  identified. 
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1.  Introduction 


Figure  1. 
Transmission 
spectrum  of 
Zn:meso-tetra  (p- 
methoxyphenyl) 
tetrabenzporphyrin 
(TBP),  taken  in  a 
1-mm  sample  with  a 
concentration  of  0.5 
g/1.  Transmission  at 
532  nm  is  71.6%. 


Programs  in  nonlinear  optical  (NLO)  materials  generally  aim  at  develop¬ 
ing  materials  that  exhibit  good  optical  quality,  high  damage  thresholds, 
large  nonlinear  behavior,  and  fast  response  times,  and  that  lend  them¬ 
selves  to  easy,  low-cost  manufacture.  Many  material  systems  have  been  in¬ 
vestigated,  but  organic  molecules  offer  perhaps  the  greatest  challenge.  In 
addition  to  the  sheer  size  of  this  group  of  materials,  many  of  the  nonlinear 
processes  that  occur  with  these  materials  are  not  well  understood.  Despite 
these  obstacles,  organic  molecules  are  attractive  for  a  number  of  reasons. 
They  are  potentially  low  cost  and  easy  to  S5mthesize.  They  can  often  be  in¬ 
corporated  into  polymers  or  tailor-made  for  specific  device  designs.  Fur¬ 
ther,  organics  have  exhibited  some  of  the  largest  nonlinear  optical  re¬ 
sponses  yet  found. 

An  important  class  of  organic  NLO  materials  is  that  of  organometallic  ma¬ 
terials,  such  as  the  phthalocyanine  and  porphyrin  materials.  These  materi¬ 
als,  which  contain  metals  surrounded  by  highly  conjugated  ji-bonds,  can 
exhibit  large  nonlinear  responses.  Previous  investigators  have  noted  that 
the  response  of  an  organometallic  material  to  wavelength,  as  well  as  its 
nonlinear  magnitude,  can  be  modified  by  variation  in  the  metal  atom  that 
resides  in  the  center  of  the  molecules  [1]  or  in  the  constituents  attached  to 
the  macrocyclic  ring.  In  this  report,  we  highlight  the  nonlinear  optical  be¬ 
havior  of  porphyrin-type  organometallic  materials,  specifically,  Zn:meso- 
tetra(p-methox5q)henyl)  tetrabenzporph5nrm  (TBP)  (fig.  1). 
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2.  Previous  Results 


Earlier  studies  on  TBP 


reported  a  second  order  hyper-polarizability,  |(y)5|/ 


of  4.8  X  esu  at  532  nm  with  picosecond  pulses  [2].  Other  researchers 
have  shown  that  the  nonlinear  transmission  in  TBP  is  fluence  dependent 
and  that  the  dominant  nonlinear  mechanism  is  excited-state  absorption  [3]. 
The  same  mechanism  has  been  found  in  other  porphyrin-type  molecular 
systems  [4]. 


In  our  previous  investigations  at  532-nm  wavelength  and  ~10-ns  pulse  du¬ 
ration  [5],  we  foimd  that  TBP  has  one  of  the  largest  reported  excited-state 
absorption  cross  sections  (<7^^  ~  240  x  10"^®  cm^)  and  excited-state  to 
ground-state  cross-section  ratios  (~28),  and  that  the  index  change  at 
532-nm  is  due  to  the  contribution  of  thermal  effects  and  an  excited-state 
population.  Because  the  excited-state  refractive  component,  (T^,  was  fovmd 
to  be  negative,  its  effect  was  added  to  the  thermal  density  change  (which  is 
also  negative);  this  result  is  unlike  previously  reported  results  in  macrocy- 
clic  organometallics.  We  have  also  shown  [6]  that  TBP  performs  well  when 
compared  to  similar  organometallic  materials  and  carbon  black  suspen¬ 
sions  (CBS)  in  an//5  optical  limiting  configuration.  The  TBP  for  the  earlier 
study  [6]  and  the  work  reported  here  was  dissolved  in  tetrahydrofuran 
(THF)  at  a  concentration  of  C  =  0.5  g/1  (where  a|532  cm“^). 

In  our  previous  work,  we  measured  the  limiting  behavior  of  TBP  with  the 
sample  placed  at  focus  in  an//  64  optical  system;  this  location  corresponds 
to  a  spot  size  (half  width  e~‘^  maximum)  at  focus,  of  21.4  pm.  The  trans¬ 
mitted  energy  versus  the  input  energy  was  measured  through  an 
apertured  detector  set  to  collect  only  the  1  /e^  cone  of  the  output  Gaussian 
spatial  profile  energy  (86.5  percent).  We  used  parameters  determined  from 
a  Z-scan  analysis  to  model  the  limiter  performance  of  this  material.  The 
limiting  data  (points)  and  the  theoretical  fits  (curves)  are  shown  in  figure  2. 
The  three  curves  indicate  linear  behavior  (dashed  line),  the  theoretical  fit 


Figure  2.  //64  limiter 
performance. 
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assuming  that  all  the  absorbed  energy  was  converted  into  heat  (solid  line), 
and  a  theoretical  fit  assuming  that  only  the  energy  from  the  linear  absorp¬ 
tion  was  converted  into  heat  (long  dashed  line).  The  data  were  found  to  fit 
the  latter  model  up  to  an  input  energy  of  ~3  |ij,  at  which  point  the  analysis 
broke  down.  As  can  be  seen  from  the  theoretical  fits,  the  best  limiting  oc¬ 
curs  when  all  the  absorbed  energy  can  be  converted  into  heat.  The  data 
did  not  fit  the  best-case  scenario  (solid  line),  possibly  because  of  other  con¬ 
tributions  not  accounted  for  in  this  model,  such  as  fluorescence,  higher 
level  transition  states,  material  impurities,  breaking  of  molecular  bonds, 
saturation,  boiling,  or  the  acoustic  wave  not  fuUy  propagating  across  the 
beam  during  the  length  of  the  pulse.  We  explore  several  of  these  scenarios 
in  this  report  and  determine  the  most  likely  explanation. 

3.  Current  Experiments 

Although  two  effects  are  causing  nonlinear  transmission  (nonlinear  ab¬ 
sorption  and  self-defocusing),  the  theoretical  heating  solid  curve  in  figure 
2  is  all  self-defocusing.  Energy  is  not  reaching  the  detector  because  a  por¬ 
tion  of  it  is  spread  over  such  a  spot  size  that  it  cannot  fit  through  the  aper¬ 
ture  about  the/-cone,  as  it  could  at  low  input  energies.  Heating  causes  self- 
defocusing  by  changing  the  index  of  refraction.  The  incident  radiation  is 
absorbed  by  the  TBP,  and  this  energy  is  returned  as  heat,  which  increases 
the  temperature  of  the  surrounding  THE  solvent.  The  heated  solvent  then 
changes  density,  which  in  turn  changes  the  index  of  refraction.  Since  the 
density  is  decreased,  the  index  is  decreased,  so  that  the  material  forms  a 
highly  aberrated  negative  lens. 

We  confirmed  the  existence  of  fluorescence  in  Zn:TBP  experimentally.  A 
time-averaged  photoluminescence  spectrum  was  measured,  which  con¬ 
firmed  that  upon  excitation  at  532  nm,  a  fluorescence  exists  at  658  nm  and, 
to  a  lesser  extent,  at  630  nm.*  Fluorescence,  which  allows  the  excited  mol¬ 
ecule  to  de-excite  without  the  release  of  heat,  could  contribute  to  the  devia¬ 
tion  of  our  data  from  the  theoretical  fit;  however,  the  observed  fluorescence 
does  not  appear  to  be  strong  enough  to  account  for  the  large  deviation  in 
figure  2. 

The  density  change  is  not  instantaneous.  It  occurs  on  the  order  of  the  speed 
of  sound.  What  is  more,  it  launches  a  sotmd  wave  that  propagates  beyond 
the  heated  region,  so  it  is  a  nonlocal  event.  This  additional  complication  to 
the  analysis  is  not  incorporated  into  the  theory  used  to  derive  figure  2.  The 
measured  spot  size  for  the  limiting  performed  in  figure  2  was  21.4  pm.  As¬ 
suming  an  acoustic  velocity  t)^ically  found  in  organic  solvents  of  v  ~  1.5  x 
10^  cm/s,  the  acoustic  wave  generated  by  the  change  of  density  caused  by 
the  temperature  change  will  spread  to  only  about  half  the  spot  size  in  9  ns. 
Therefore,  the  heat-generated  density  change  in  our  model  is  likely  an 
overestimate,  as  the  experiment  was  conducted  in  a  transient  regime.  Since 
orvly  the  centralmost  portion  of  the  beam  was  affected  by  the  thermally  in¬ 
duced  index  change,  this  arrangement  reduced  the  refractive  optical  limit¬ 
ing  ability  of  TBP. 


*  Private  conversation  with  Mary  Tobin,  ARL,  Adelphi,  MD. 
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We  devised  a  new  experimental  arrangement  for  the  limiting  experiment 
with  an  // 14  optical  input  {w^  =  4.96  fi.m  at  focus)  to  allow  the  acoustic 
wave  the  time  to  traverse  the  entire  laser  beam  during  the  pulse.  In  this  ar¬ 
rangement,  the  entire  thermal  response  created  by  the  incident  laser  pulse 
should  enhance  the  limiting  results,  and  the  density  change  should  reach 
steady  state.  We  used  a  lOO-pm  sample  length  to  satisfy  the  "thin"  cell  cri¬ 
terion,  L  «  Zp,  that  is  required  for  the  modeling  we  employ.  In  our  case,  the 
Rayleigh  range,  z^  =  tzw^/X,  was  141  |im,  where  zu^  =  spot  size  at  focus.  In 
this  optical  setup,  the  acoustic  wave  could  fully  traverse  the  laser  pulse, 
but  at  intensity  levels  much  higher  than  in  the  previous  experiments. 

We  performed  a  Z-scan  analysis  under  the  new  focusing  conditions  first  to 
obtain  the  appropriate  parameters  for  the  limiting  model  [7].  Our  work 
with  this  material  used  a  Continuum  Q-switched,  injection-seeded 
Nd:YAG  laser  operating  at  10  Hz.  The  laser  was  frequency  doubled  to  a 
wavelength  of  532  nm  and  operated  in  the  fundamental  transverse  electro¬ 
magnetic  (TEMgg)  mode.  We  used  a  dual-leg  Z-scan  setup  to  reduce  noise 
in  the  data.  As  stated  above,  the  beam  radius  zo^  was  4.96  |im  at  focus.  The 
temporal  pulse  width  at  full  width  at  e~^  intensity  was  =  13  ns.  To  deter¬ 
mine  the  excited-state  cross  section  under  our  experimental  conditions,  we 
assumed  a  three-level  model  (fig.  3),  with  (1)  optical  transitions  from  the 
ground  state  to  the  excited  singlet  state  Sj,  (2)  intersystem  crossing  to  a 
triplet  state  T^,  and  (3)  optical  transitions  from  T^  to  an  upper  level.  We 
also  assumed  that  no  measurable  saturation  or  diffusion  effects  occurred 
during  the  pulse,  and  that  the  intersystem  crossing  rate  was  fast  compared 
to  the  pulse  duration  «  t  ).  Equation  (1)  describes  the  change  in  in¬ 
tensity  with  distance  through  me  material: 

f  =  ,  (1) 

where  is  the  number  density  of  ground  state  charges,  Nj-^  is  the  num¬ 
ber  density  of  charges  in  the  lowest  triplet  state,  the  cr's  are  the  cross  sec¬ 
tions  for  absorption  in  the  ground  and  first  triplet  state,  and  I  is  the  beam 
intensity.  Note  that  this  equation  assumes  that  the  spot  size  does  not 
change  within  the  sample;  i.e.,  diffraction  is  small  and  ignored.  Solving 


Figure  3.  Three-level 
model  used  to 
describe  Zn:TBP. 
Diagram  on  right 
represents  a  Zn:TBP 
molecule.  Cross 
sections  and 
relaxation  times  are 
given.  K's  are  decay 
rates.  Subscripts: 
isc  =  intersystem 
crossing;  es  =  excited 
state;  gr  =  ground,  Tl, 
SI,  Tl,  etc,  =  singlet, 
triplet  states. 


Singlet  states 


Triplet  states 


Absorption  cross  section,  a  (cm^) 
8x10-1® 

a,,  240  X  10-1® 


Relaxation  time,  t 
x,j,,  «10ns 

»10ns 
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this  equation  for  the  fluence  and  energy  density,  and  integrating  over  the 
spatial  extent  of  the  beam,  we  may  write  the  normalized  energy  transmis¬ 
sion  T  through  a  material  of  length  L  as 


% 

l+x2 


where 


'^0"  2im  ^  ' 

X  =  z/ Zg,  z  is  the  position  of  the  sample  along  the  optical  path  (z  =  0  at  fo¬ 
cus),  Oq  =  =  0),  and  Fg  is  the  input  fluence  as  a  function  of  r  and  po¬ 

sition  z.  The  above  equation,  which  assumes  an  infinite  number  of  charges 
available  in  the  groimd  state,  can  be  used  to  fit  the  open-apertured  Z-scan 
data  and  get  the  value  for  ^g.  From  ^g,  can  be  obtained,  where  O  is  the 
triplet  yield  and  can  vary  between  0  (no  excited  charges  transfer  from  the 
Sj  to  Tj  state)  and  1  (all  the  excited  charges  transfer  from  the  Sj  to  state). 
A  typical  open-aperture  (where  all  the  transmitted  light  is  collected)  Z- 
scan  curve  and  fit  are  shown  in  figure  4. 

The  apertured  Z-scan  data  contain  both  the  nonlinear  absorption  and  re¬ 
fraction  components  of  the  material's  nonlinear  behavior.  Using  an  aper¬ 
ture  that  transmits  only  25  percent  of  the  beam  (S  =  0.25)  in  the  far  field  (a 
distance  d  from  the  sample,  where  d  »  Zg),  one  could  iteratively  fit  the  nor¬ 
malized  energy  transmission  to  determine  the  excited-state  population  re¬ 
fractive  index  cross  section,  <7^,  assuming  a  value  for  the  heat-induced  in¬ 
dex  change.  Only  two  contributions  were  assumed  to  add  to  the  nonlinear 
refraction — one  derived  from  a  change  in  density  due  to  heating  of  the 
sample,  and  the  other  from  the  excited-state  population  contribution  to  the 
index. 


Figure  4.  Tjrpical  open- 
apertured  Z-scan. 

Solid  circles  represent 
data  and  solid  line 
represents  a  fit  using 
eq  (2). 
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Figure  5.  Open- 
apertured  Z-scan  for  a 
saturated  condition, 
energy  £  =  69  nj.  Solid 
circles  represent  data 
and  solid  line 
represents  a  fit  using 
eq  (2). 


4.  Results  and  Discussion 

The  open-aperture  Z-scan  data  from  the// 14  optical  setup  was  fit  by  equa¬ 
tion  (2).  At  energies  above  7  nJ,  this  equation  could  not  exactly  fit  the  data 
(see  fig.  5).  Experimentally,  the  Z-scan  curves  were  wider  than  expected — 
i.e.,  equation  (2)  predicts  a  narrower  profile  based  on  the  measured  spot 
size  of  the  optical  system.  By  varying  in  equation  (2)  to  achieve  its  best 
fit  to  the  data,  we  found  that  the  predicted  was  fluence  dependent.  As 
the  fluence  increased,  decreased.  At  input  energies  below  ~1  nJ  (2.6 
mj/ cm^),  remained  constant  at  the  earlier  reported  value  (~240  x  10“^® 
cm^)  but  with  a  larger  imcertamty.  The  excited-state  cross  section  should 
be  a  constant;  therefore  its  fluence  dependence  indicates  a  breakdown  in 
the  simple  model  assumed  here,  at  least  above  ~1  nJ.  The  most  obvious 
omission  was  not  incorporating  saturation  of  the  ground  state  into  the 
model.  If  this  saturation  is  incorporated  for  the  fluence,  equation  (2) 
becomes 


where  Fg  =  saturation  fluence,  defined  as 


F  -  h(o 

s  (bcr 
^  gr 


(4) 


Unlike  equation  (2),  this  equation  cannot  be  solved  in  closed  form  and 
must  be  numerically  integrated. 


While  Fg  could  not  be  determined  exactly  because  of  the  uncertainty  in  the 
triplet  yield,  we  were  confident  that  the  fluences  encoimtered  during  the 
optical  limiting  experiment  were  well  above  Fg,  which  is  on  the  order  of 
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0.05  to  0.15  J/cm^.  While  adding  saturation  of  the  ground  state  allowed  for 
a  more  accurate  fit,  perfect  fits  to  the  theory  could  not  be  obtained.  One 
possible  explanation  for  this  could  be  triplet-state  saturation. 

Once  again,  we  performed  optical  limiting  as  described  above,  only  using 
the  new// 14  optical  arrangement.  The  data  and  theoretical  fits  are  shown 
in  figure  6.  We  allowed  the  heat  to  fully  cross  the  beam  during  the  pulse 
over  the  fluence  range  that  we  used  in  our  previous  work;  however,  the 
limiting  curves  did  not  significantly  improve.  Changing  the  spot  size  did 
not  result  in  the  expected  thermal  contribution,  but  this  could  be  explained 
by  saturation  of  the  groimd  state.  Saturation  prevents  the  data  from  limit¬ 
ing  as  well  with  increasing  input  energy  for  two  reasons.  First,  the  absorp¬ 
tion  does  not  grow  as  large  as  expected;  second,  the  smaller  amount  of 
absorption  means  that  the  solvent  is  not  heated  as  much,  so  that  the  self- 
defocusing  is  not  as  large.  It  is  interesting  to  note  in  figure  6  that  limiting 
above  0.2  pj  is  better  than  would  be  predicted  based  on  nonlinear  absorp¬ 
tion  and  nonlinear  refraction,  with  heat  coming  only  from  linear  absorp¬ 
tion.  Since  saturation  of  the  absorption  most  likely  is  occurring,  the 
absorption  does  not  increase  as  predicted  here,  but  some  heating  is 
occurring  from  the  excited  state  as  well — ^just  not  as  much  as  expected 
based  on  the  xmsaturated  case.  We  have  not  incorporated  saturation 
effects  into  the  Umiting  model  to  check  this  h5^othesis  but  expect  to  in  the 
near  future. 


Figure  6.//14  limiter 
performance. 
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5.  Conclusions 


To  accurately  model  this  molecule,  one  cannot  ignore  saturation  effects.  In 
fact,  our  data  indicate  that  not  only  is  groimd-state  saturation  significant, 
but  saturation  of  the  triplet  states  may  also  be  occurring.  Work  is  continu¬ 
ing  so  that  we  can  better  vmderstand  this  behavior  and  quantify  these 
saturations. 
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